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Abstract 

We  consider  the  problem  of  estimating  the  rate  of  a  doubly-stochastic,  time-space  Pois¬ 
son  process  when  the  observations  are  restricted  to  a  region  D  C  Ft2.  In  the  general  case,  we 
obtain  a  representation  of  the  minimum  mean-square-error  (MMSE)  estimate  in  terms  of  the 
conditional  characteristic  function  of  an  underlying  state  process.  In  the  case  D  —  IR2,  we 
extend  a  known  result  to  compute  the  MMSE  estimate  explicitly.  For  a  special  form  of  the 
rate  process,  a  well-defined  integral  equation  is  presented  which  defines  the  linear  MMSE  esti¬ 
mate  of  the  rate. 


Key  Words-,  doubly-stochastic,  time-space  Poisson  process,  MMSE  estimate,  linear  MMSE 
estimate,  likelihood  ratio. 
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I.  Introduction 

We  consider  a  doubly-stochastic,  time-space  Poisson  process  N°  with  intensity  function 
X({,  r)  =  /  (f ,  r  -H(t)xt),  where  t  >0  and  r£  H2.  Here,  /  is  a  known,  deterministic 
function;  xt  €  IR”  is  the  solution  of  an  Ito  stochastic  differential  equation,  and  H(t)  is  a 
known,  deterministic,  1R2Xb -valued  function.  The  process  N°  under  consideration  counts 
events  which  occur  in  all  of  1R2;  however,  suppose  that  only  those  events  which  occur  within  a 
region  D  C  ]R2  can  be  observed.  We  wish  to  compute  minimum  mean-square-error  (MMSE) 
estimates  of  \(t  ,  r  ),  given  our  limited  observations.  In  the  general  case,  D  1R2,  we  obtain 
a  representation  of  these  estimates  in  terms  of  the  conditional  characteristic  function  of  xt  . 

-  — r ' R  (t  r'r 

When  D  =  IR2,  and  /  (f ,  r  )  =  e  2  .for  some  deterministic  matrix  /?((),  we  extend 

a  result  of  Rhodes  and  Snyder  [l]  to  compute  the  MMSE  estimate  of  X(f ,  r )  explicitly.  We 
also  consider  linear  estimates  of  X(f  ,  r)  for  the  same  choice  of  /  when  D  )R2.  These 
filtering  problems  are  frequently  encountered  in  optical  communication  systems  [2,  3],  particu¬ 
larly  in  the  context  of  hypothesis-testing;  this  issue  is  discussed  in  Section  V. 

H.  Probabilistic  Setting 

Let  J32  denote  the  Borel  subsets  of  IR2.  Next,  if  I  is  any  interval  of  IR,  let  IB  (I)  denote 
the  Borel  subsets  of  /.  We  define  IB  (I  )(£>  IB2  to  be  the  smallest  a-field  containing  all  sets  of 
the  form  E  X  A  ,  such  that  E  E  IB  (I)  and  A  €  IB2.  Let  (  B  ,  F  ,  P  )  be  a  probability  space 
on  which  we  let 

N°  =  {  N(B)  :  BE  B(0.oo)®B2  }  , 

be  a  time-space  point  process.  Sometimes,  N°  is  called  a  random  point  field  or  a  random 
measure.  Here,  this  means  that  with  each  B  €  IB  (0,oo)®  IB2,  we  associate  a  nonnegative,. 
integer-valued  random  variable,  N  (B  )  —  N  (w,  B  );  in  addition,  for  each  w  (E  O,  N  (  w,  •  )  is 
assumed  to  be  an  integer-valued  measure  on  IB  (0,oo)®  IB  V  We  let  Ft  represent  the  times 
and  locations  at  which  points  have  occurred  up  to  and  including  time  t  .  More  precisely,  let 
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F0  denote  the  trivial  cr-fleld,  and  for  t  >0,  set 

Ft  =  a{  N(B)  :  Be  B(0,t}®B*  }  . 

Now,  let  D  be  a  Borel  subset  of  1R2.  We  take  G0  to  be  the  trivial  er-field,  and  for  t  >0,  we 
set 

Gt  =  er{  N(  Bp)  {  (0,oo)  X  D  }  )  :  Be  B(0,t}®B'2  }  . 

Note  that  €  t  represents  the  history  of  the  point  process  restricted  to  the  region  D  ,  up  to 
time  t  .  We  shall  refer  to  G  t  as  our  “observations  up  to  time  t On  the  same  probability 
space,  (  f2  ,  F ,  P  ),  let  Ar  be  an  n -dimensional  Gaussian  random  vector  with  known  mean, 
m  ,  and  known,  positive-definite  covariance,  S  .  Let  {  vt  ,  t  >0  }  be  a  standard  Wiener  pro¬ 
cess  independent  of  X .  We  let  the  n  -dimensional  process  {  xt ,  t  >0  }  be  the  solution  to  the 
Ito  stochastic  differential  equation 

dxt—F(t)xtdt  +  V(t)dvt  ;  x0  —  X  .  (1) 

Here  F  and  V  are  known  matrices  with  appropriate  dimensions.  We  also  assume  that  F  and 
V  are  piecewise-continuous  so  that  a  unique  solution  of  (l)  exists  (see  Davis  [4],  pp.  108-111). 
Let 

JT0  =  g{  x,  ,  0<s  <oo  }  . 

For  t  >0,  let  denote  the  smallest  cr-fleld  containing  Ft  JT0.  We  write  this  symboli¬ 
cally  as 

X,  =  F,  v  X0  ;  t  >0. 

We  shall  assume  that  N°  is  an  {Xt  }-doubly-stochastic,  time-space  Poisson  process,  with  X0- 
measurable  intensity  (see  Bremaud  [5],  pp.  21-23  and  233-238) 

X(f ,  r  )  =  /  (f ,  r  -  H (f  )x, )  , 

where  t  6  (0,oo),  r  £  IR2,  and  xt  is  defined  by  (1).  Assume  that  FI :  (0,oo)  — *  fft2Xn  and 
/  :  (O.oo)XR2  — ^  (O.oo)  are  deterministic  and  known.  We  further  assume  that  the  function 
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/*(* )  =  /r2  f  (t ,  r  )  dr  (2) 

is  finite  for  all  t  <oo.  This  means  that  for  each  f  >0,  the  process 

N‘  1  {  N(B)  :  Be  B(t ,  oo)®B2  } 

is  a  Poisson  random  field  under  the  measure  P(  •  |  Xt ),  with  rate  X(s  ,  r ),  where  e  £  (t ,  oo), 
and  r£E2.  This  implies  the  following.  First,  for  B  E  B  {0,oo)®  B2,  let 
MB)  t  L  \(s  ,  r)  dr  ds  ;  then  if  B  £  B  (t ,  oo)®  B2  and  n  is  an  arbitrary,  nonnegative 
integer, 

P(  N(B)  =  n  |  J()  =  MB  ^  e  -A(B)  ,  (3) 

n  ! 

and  hence,  for  6  E  R, 

E  [  I  Xt  ]  =  expl(e"-l)A(B)l  .  (4) 

The  second  implication  is  that  if  and  S2  are  disjoint  sets  in  JB(f ,  oo )®JB2,  then  the  ran¬ 
dom  variables  N  (  B  t  )  and  N(  B2)  are  independent  under  the  measure  P  (  •  |  JT, ). 

Notation.  We  let  N0  =  0  and  for  t  >0,  Nt  =  N(  (0,t]XD  ). 

III.  Nonlinear  Filtering  Results 

We  first  establish  some  notation  in  order  to  state  our  results  more  compactly.  Let 
Pt(x),  i£E’,  denote  the  (regular)  conditional  probability  of  xt  given  Gt  .  Let  t/>,  (r/), 
rj  E  IP"  ,  denote  the  conditional  characteristic  function  of  xt  given  € ,  : 

(V)  =  E  [  e  |  <D  t  ]  =  /  e  dPt  {x  )  ■  r?  £  Rn. 

Next,  let 

X(t ,  r)  t  E  [  X(f  ,  r )  |  Gt  ]  =  E  [  /  (f  ,  r  -  H (t  )*,  )  |  C,  ]  , 


and 
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)(t,6)  t  f  2\(t,  r)  e>',,r  dr  ;  9  E  1R2. 

We  also  set 

F(t,6)  t  f  f  (t,r)e>’lr  dr  . 
Theorem  1.  Under  the  foregoing  assumptions, 

'l(t,9)  =  F(t,9)^l(H(t)'6)  . 

Proof.  Observe  that. 

=  f  „E  [  f(t,  r  -H(t)Z')  j  Ct  }eie’r  dr 

R“ 

=  fn2  JK„  /  (t  ,  r  -  H  (t  )x  )  dPt  (x  )  e  dr  . 

By  Fubini's  Theorem, 

~l(t,e)  =  f  f  f  ( t,r  -H(t)x)  e  **>T  dr  dPt{x) 

R  R 

=  F(t,  6  )  f  dPt(x) 

=  F(t,8)J  «/(»(*)'#)’»  ip(x) 

■'r" 

=  F(t,d)rft(  H(t)'6). 


Theorem  2.  i/  D  =  IR2,  and  if 


f  (t  .  r)  —  e 


--r’tUty 


QED 


(5) 


for  some  deterministic,  positive- definite  matrix  R  ( t ),  then 
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X(t.  r)  t  E  [  X(f ,  r  )  |  Ct  } 

=  E  [  f  (t,  r  -H(t)xt)  j  Ct] 

=  exp[-j-(r-H(t)x,Y  Qf1  ( r-H(t)x  ,)  J  , 

vdet  Qt  2 

where 

xt  ~  E  [  xt  |  <?f  ]  , 

E(  =  E  [  (x,  -  x  ,  )(x(  -  x  ( )'  |  G ,  ]  >  0  ,  P  -  a.s.  , 

Qt  t  H(t)btH(t)'  +  R(t). 

and 

dx  t  =  F(t)r  t  dt  (6) 

-f  /r8E, _//(<-)'  g«:‘  (r-mt-)z,_)N(dtX  dr);  x0  =  m  , 


rfSt  =F(OS,  dt  +  E(F(t)'  dt  +  V(t)V(t)'  dt  (7) 

-  E,  _//(<-)'  Qd1  H(t-)Zt_N(dt  X  1R2  )  ;  E0  =  5. 

Proo/.  First,  since  P  =  1R2,  G,  =  F, .  Next,  in  [1]  it  is  proved  that  the  conditional  density 
of  x,  given  Ft  is  Gaussian  with  conditional  mean  x  (  and  conditional  covariance  E,  (which  is 
positive  definite  almost  surely  because  of  the  assumption  that  S  is  positive  definite)  satisfying 
(6)  and  (7)  above.  So, 


Next,  from  equation  (5),  it  follows  that 


F  (t  ,  9  )  —  2n  Vdet  R  (t )  e  2 

Hence,  by  Theorem  1, 

~l(t  .6)  =  2tt  Vdet  R(t)  - ,'Q>t  . 

Taking  inverse  Fourier  transforms,  we  see  by  inspection  that 
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When  D  ^  ]R2,  or  equation  (5)  does  not  hold,  ipt  (r?)  is,  in  general,  not  known.  This  has  led  us 
to  consider  linear  estimates  of  \(t ,  r ).  We  discuss  this  in  the  next  section. 


IV.  Linear  Filtering  Results 

We  call  \i(t  ,  r)  a  linear  estimate  of  X(<  ,  r)  given  €  t  ,  if  \i  can  be  written  in  the 

form 

\L(t  ,  r)  =  /  1  fD  h(t ,  r  ;  r,  p)  [  N  (d  r  X  d  p)  -  X(r,  p)  dr  dp)  +  h0(t ,  r  )  ,  (8) 

where,  h  and  A0  are  deterministic,  and  \{t ,  r )  =  E  [  X(f ,  r )  ].  We  wish  to  choose  h  and  h0 

to  minimize 

E  [  jX(f,  r)  -  \L(t,  r )  |  2  )  .  (9) 

Lemma  1.  (Grandell  (6]  ).  Let\L(t,  r)  be  given  by  (8).  Under  the  conditions  outlined  in 
Section  II,  the  quantity  in  (9)  will  be  minimized  if  h0(t  ,  r  )  =  \(t ,  r  ),  and  if  h  satisfice 

r{t,r-,  r,  p)  =  J  ‘  /  h  (t  ,  r  ;  a,  ?)!>,  f  ;  r.  p)  d  q  d  a  +  h  (t  ,  r  ;  r,  p)\(r,  p)  .  (10) 

where 

r(t  ,  r  ;  r,  p)  =  cov  [  X(f  ,  r  ),  X(r,  p)  \  . 

With  Lemma  1  in  mind,  we  state  our  Theorem  3. 

Theorem  3.  ///(<,  r)  is  given  by  (5),  and  the  conditions  outlined  in  Section  II  hold,  then 
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\(t,  r)  =  ^EU!Mexvl-±(r-H(t)x(t))'  Q  (l  )_1  (r-H(t  )x{t ))  ]  ,  (n) 

Vdet  Q  (t  )  2 


where 


X(t)  t  E  \xt]  , 

£(f )  =  cov  [ xt  ]  , 

<9(0  =  H(t)-z(t)H(ty  +  R(t). 


Furthermore, 


r  (t 


r,  p)  +  \(t  ,  r  )\(r,  p) 


V 


det  R  t  )  det  i?  (r) 
det  <9  (<  .  t-) 


X 


(12) 


1 

r 

H{t )  o 

x(t ) 

* 

r 

H(t )  o 

x{t) 

exp[  -  — ( 

P 

— 

0  H^t\ 

I(t) 

)'  Q(t,  r)-1  ( 

p_ 

o  H  (t) 

?  (r)  _ 

where 


S  (1  ,  r)  =  cov  (  art  ,  xr  ]  , 

and 


Q(t.r)  t 


<9(0  H(t)nt.T)H(T)' 
H(t)Z(t,  t)H(ty  Q  (r) 


Proof.  For  completeness,  we  make  the  following  observations.  Recall  that 

dxt  —F(t)xtdt  +  V(t)dvt  ;  x0  —  X.  (13) 

Let  <f>(f  2.  f  j)  be  the  transition  matrix  corresponding  to  F  (t ).  Then 

x  (t )  =  ,  0 )m  ,  (14) 


and 
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m!n(f ,  r ) 

E(f,  r)  =  *(t.  0)S4>(t,0)'  +  fg  *(t.a)  V(a)V(sy  *(r,sy  ds  . 

Note  that  E(f )  =  £(f ,  f ). 

To  compute  X(f  ,r)=E[  X(f ,  r )  ),  observe  that  xt  is  Gaussian  with  mean  x  (f )  and 
covariance  £(f ).  By  considering  the  proofs  of  Theorem  1  and  Theorem  2,  equation  (11)  is 
immediate. 


The  computation  of  (12)  is  similar,  but  requires  some  judicious  preliminary  arithmetic. 
First,  observe  that  r(<  ,  r  ;  r,  p)  +  X(f  ,  r  )X(r,  p)  is  just  another  way  of  writing 
E  [  X(t ,  r  )X(r,  p)  ].  Next,  rewrite  \{t ,  r  )X(r,  p)  as 


exp[  - 


T( 


r 

P 


H(t) 

0 

)' 

R  ( t  )-’  o 

( 

r 

//(f)  0 

0 

H(t) 

X r. 

0  /?  (r)-1 

P 

0  H(t) 

xr_ 

which  is  equal  to 


exp[  -  — ( 


M 

//(f)  o 

xt 

)' 

fi(f) 

0 

( 

H 

H  (i )  o 

xt 

k 

o  H(t) 

x  r 

0 

/?(r) 

0  ff(r) 

X  T  _ 

fx‘ 


)].  (15) 


is  a  Gaussian  random  vector  with  mean, 


Because  {  x,  ,  t  >0  }  is  a  Gaussian  process,  | 

.  By  the  same  reasoning  used  to  deduce  (11),  (12) 


*(0 
x  (r) 


and  covariance 

also  follows. 


E (t )  £(t ,  t-) 

E(r,  ( )  E(r) 


QED 


Remark.  In  equation  (10),  if  we  regard  t  and  r  as  fixed,  and  divide  through  by  X(r,  p),  then 
the  result  has  the  form  of  the  Fredholm  equation 

g  =  Bh  -7-  h  , 

for  known  function  g  ,  known, operator  B  ,  and  unknown  function  h  . 
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V.  Discussion 

The  filtering  problems  considered  above  often  arise  in  the  design  and  implementation  of 
receivers  for  optical  communication  systems.  Typically,  a  binary  message  source  is  used  by  a 
transmitter  to  select  the  modulation  of  the  intensity  of  a  laser  beam  in  accordance  with 

whether  a  “0”  or  a  “1”  is  to  be  sent.  The  laser  beam  travels  to  a  receiver  and  strikes  its  pho¬ 

todetector.  We  assume  that  the  laser  beam  has  an  intensity  profile  of  the  form 

".  (<  ){  (t  ,  r)  ;  i  =  0,  1. 

Here,  is  a  known,  deterministic  function,  where  i  =  0  or  1  has  been  selected  by  the 

transmitter. 

We  model  the  surface  of  the  receiver's  photodetector  as  IR2.  If  the  receiver,  for  exam¬ 
ple,  is  subject  to  vibrations,  the  center  of  the  spot  cf  laser  light  may  wander  randomly  over 
the  photodetector  surface  [2],  We  assume,  as  in  [2],  that  the  center  of  the  spot  of  laser  light  is 
given  by  H(t)xtE  IR2.  The  output  of  photoelectrons  from  the  photodetector  is  modeled  by 
the  process  N°,  with  stochastic  intensity  now  given  by 

X,(< .  r  )  =  Ui(t)f  (t ,  r  -  H(t  )xt )  .  (16) 

Of  course,  an  actual  photodetector  does  not  have  an  infinite  photosensitive  surface.  We 

account  for  this  fact  by  assuming  that  only  those  photoelectrons  which  occur  in  a  region 
D  C  IR2  are  observed.  For  example,  in  this  setting.  D  might  be  a  square  or  a  circle  centered 
at  the  origin.  After  observing  photoelectrons  occurring  in  D  during  some  time  interval 
[0,  T],  a  decision  as  to  whether  a  “0”  or  a  “1"  was  sent  has  to  be  made  based  on  one  of  the 
estimates  X,(f,  r)  or  X,-  i(f,  r  ).  As  an  example  of  a  decoding  scheme,  we  could  use  the 
likelihood  ratio  test 

> 


l, 
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to  make  the  decision,  using  the  minimum  probability  of  error  cost  criterion  and  assuming 
equiprobable  hypotheses  (see  Snyder  [3],  section  2.5).  The  likelihood  ratio,  Lr ,  is  given  by 
(see  Snyder  [3],  pp.  471-476) 


nT  .  T 

n  Xt(fy  ,  rj  )  exp[  -  fo  fD  X,(.s  ,  r  )  dr  ds] 

Lt  —  JLri _ 

„  T 

n  Xo(«j  .  rj )  exp[  -  Jo  fD  x0(s  ,  r  )  dr  ds) 


(17) 


where  t:  and  r;-  are  respectively  the  time  and  the  location  of  the  j  th  photoevent  in  the 
region  D  ,  and  we  adopt  the  convention  that  when  NT=  0,  the  factors  proceeding  exp  in 
equation  (17)  are  taken  to  be  unity.  Here,  of  course, 

X,(f,  r)  *  E  [X,-(t,  r)  |  Gt  ]  ;  t  =0,  1. 

Now,  using  (16),  (17)  simplifies  to 

Nt  "i(0  .7* 

Lt  =  n  -LJ-  exp[  -  J  J  X,(«  ,  r  )  -  X0(«  ,  r  )  dr  ds  ]  .  (19) 

}  =1  ) 


In  the  general  case,  D  ft2,  X,(f ,  r)  is  not  known,  and  hence,  L?  cannot  be  com¬ 
puted.  However,  when  D  —  1R2,  it  turns  out  that  we  do  not  need  to  know  X,-  (( ,  r )  in  order 

to  compute  Lr  .  Observe  that  if  D  =  IR2,  then  : 

JD  ~xi(*  >  r  )  -  X0(s  ,  r  )  dr  —  E  [  J^2  X,(«  ,  r  )  -  X0(s  ,  r  )  dr  |  €,  ] 

=  E  [  (  i/,(s  )  -  t/0(s  ))  J  „  f  (a  ,  r-H(s)X')  dr  \  € ,  ] 

—  E  [  (  p,(s  )  -  i/0(s  )  )  n(s  )  ]  C,  j  (20) 

=  /<(■“>'  )  [  ^[(s  )  -  "o(*  )  1- 


In  equation  (20)  we  used  the  fact  that  for  all  r0  E  R2, 

)  =  J  „  f  (s  ,  r)  dr  =  J  f  (s  ,  r  ~  r0)  dr. 
Thus,  when  D  =  IR2,  (19)  becomes 
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=  n  — 77  T  exp[  -  L  f*(s  )t  *T(«  )  -  ^o(«  )  ]  *  ]  •  (21) 

;=1  u0\lj  ) 

With  (21)  in  mind,  consider  the  following  theorem. 

Theorem  4.  The  random  field 

M‘  *  {  N(  Ex  R2  )  :  Ee  B(t,  oo)  }  , 

is  independent  of  the  cr-field  JT,  . 

Proof.  To  prove  that  M‘  is  independent  of  JT(  ,  it  is  sufficient  to  show  that  the  conditicn;il 
characteristic  function  of  N(  E  X  IR2  )  is  deterministic  for  EE  B  (t ,  oo).  Now,  it  follows 
immediately  from  the  assumption  that  N°  is  an  {■*,}  -doubly-stochastic,  time-space  Poisson 
process,  that  for  9  E  R, 

E  [  e  JSA'(  £  xr2)  j  X,  ]  =  exp[  (e  "  -  1)  f£  X,-  (s  ,  r  )  dr  ds  } 

=  exP[  1)  f£  )  /r2  /  (*  .  r-H(s  )x,  )  dr  ds  ] 

=  exp[  (ejt-  1)  f£  i/i(s  )  //(s  )  ds  ]  . 

Hence  Mf  is  independent  of  X,  . 

QED 

It  follows  from  equation  (21)  and  Theorem  4  that  for  all  t  >0,  the  random  variable  Lt  is 
independent  of  the  cr-field  X(  . 

If  we  replace  equation  (l)  by 

dxt  —F{t)xtdt  +  G(t)utdt  +  V (t  )dvt  ;  x0  —  X,  (22) 

where  {  ut ,  (  >0  }  is  predictable  with  respect  to  {  C,  ,  (  >0  }  and  G  (t )  Is  a  known  matrix 

with  appropriate  dimensions,  then  most  of  the  above  results  hold  with  only  minor 
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modifications.  The  term  G  (f  )u(  in  (22)  is  interpreted  as  a  control  signal  driven  by  the  outr 
put  of  the  photodetector.  Since  H(t)xt  represents  the  center  of  the  spot  of  laser  light  strik¬ 
ing  the  receiver,  one  might  try  to  use  G  (t)ut  to  drive  xt  to  the  origin.  This  problem  is 
addressed  in  [l].  If  (l)  is  replaced  by  (22),  Theorem  1  is  unchanged.  Theorem  2  still  holds 
except  that  equation  (6)  must  be  replaced  by 

dx  (  =  F  (t  )x  i  dt  +  G(t)utdt 

+  f  St_H(t-y  QtZ1  (r-H(t-)x  ,_)  N  (dt  X  dr)  ;  x0=m  . 

Lemma  1  is  unchanged,  and  if  ut  —  u  (t )  for  some  deterministic  control  {  v  (t ),  t  >0  },  then 
Theorem  3  holds;  of  course,  (13)  becomes  (22)  and  (14)  is  replaced  by 

f 

z(t)  —  <P(t,0)m  +  Jq  <t>(t  ,  s  )G  (s  )u  (s  )  ds  . 


In  addition,  the  results  of  the  preceding  paragraphs  of  Section  V,  including  Theorem  4,  are 
unchanged  by  substituting  equation  (22)  for  equation  (1).  Note  also  that  since  Gt  C  JCj ,  and 
Lt  is  independent  of  Xt  when  D  =  1R2,  it  follows  that  LT  is  independent  of  the  control  law 
{«t,0<(<T  }  when  D—  1R2.  This  implies  that  the  probability  of  a  decoding  error 
corresponding  to  the  likelihood  ratio  test  preceding  equation  (17)  is  not  a  function  of  the  con¬ 
trol  law  {  ut,0<t<T  }  when  D  —  R2.  In  this  sense,  all  controls:  are  optimal,  when 
D  =  R2.  In  general,  when  D  7^  1R2,  this  is  not  to  be  expected. 
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